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Abstract 

Using a calibration method, we prove that, if w is a function which satisfies 
all Euler conditions for the Mumford-Shah functional on a two-dimensional 
open set Q, and the discontinuity set of w is a segment connecting two 
boundary points, then for every point (xo,ya) of Q there exists a neigh- 
bourhood U of (xo,yo) such that w is a minimizer of the Mumford-Shah 
functional on U with respect to its own boundary values on dU . 
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1 Introduction 



This paper deals with Dirichlet minimizers of the Mumford-Shah functional (see || and |^|) 

^uix^dxdy + UHSu), (1.1) 



where fi is a bounded open subset of R with a Lipschitz boundary, H. 1 is the one-dimensional Hausdorff 
measure, S u is the set of essential discontinuity points of the unknown function u , while Vw denotes its 
approximate gradient (see |2| or 



A Dirichlet minimizer of ( |l.l| ) in £1 is a function w which belongs to the space SBV(fl) of special 
functions of bounded variation in (see or Q) and satisfies the inequality 

\7w(x,y)\ 2 dxdy + H 1 {S W ) < [ \\7u{x,y)\ 2 dxdy + H 1 {S U ) 
n Jn 

for every function u S SBV(fl) with the same trace as w on dQ. 



Suppose that w is a Dirichlet minimizer of (1.1) in £1 and that S w is a regular curve. Then the 
following equilibrium conditions are satisfied (see |5[] and ||): 

(a) w is harmonic on Q \ S w ; 

(b) the normal derivative of w vanishes on both sides of S w ; 

(c) the curvature of S w is equal to the difference of the squares of the tangential derivatives of w on 

both sides of S w . 

Elementary examples show that conditions (a), (b), and (c) are not sufficient for the Dirichlet minimality 
of w. 

In this paper we prove that, if S w is a straight line segment connecting two points of dCl, and the 
tangential derivatives d T w and d 2 w of w do not vanish on both sides of S w , then (a), (b), and (c) imply 
that every point (xo,yo) in Q has an open neighbourhood U such that w is a Dirichlet minimizer of 



( |L1 ) in U . In other words, under our assumptions, conditions (a), (b), and (c) are also sufficient for the 
Dirichlet minimality in small domains. We hope that our proof will be useful in the future to achieve the 
same result without our special assumptions on S w . 

The proof is obtained by using the calibration method adapted in (!]] to the functional (1.1). We 



construct an explicit calibration for w in the cylinder C/xR, where U is a suitable neighbourhood of 
(xo,yo). This construction is elementary when (xo,yo) ^ S w (sec |lj]), so we consider only the case 
(xo,Vo) € S w . 

The plan of the paper is the following. In Section 2 we fix the notations and we recall the main result 
of In Theorem 3.1 we consider the special case of the function 



x if y > 0, 
—x if y < 0, 



w(x,y) := 

and give in full detail the first example of a calibration for a discontinuous function which is not locally 



constant. In Theorem 3.2 we adapt the same construction to the function 

w(x,y) 



x + 1 if y > 0, 
x if y < 0. 
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In Section 4 we consider the general case of a function w satisfying (a), (b), and (c) and with S u 
{(x,y) G SI : y — 0} . If S w is connected, only two situations are possible: 



d x w{x, 0+) = — d x w(x, 0— ) on S w , 
d x w(x, 0+) = d x w{x, 0— ) on S w . 



(1.2) 
(1.3) 



The former case (1.2) is studied in Theorem |4.l| by a suitable change of variables and by adding two new 
parameters to the construction used in Theorem [O. The minor changes for the case (1.3) are considered 



in Theorem 4.2 



2 Preliminary results 

Let SI be a bounded open subset of R 2 with a Lipschitz boundary and let 

Sl a = {(x,y) e Q : y £ 0}, S = {(x,y) e SI : y = 0}. 
For every vector field : St x R — ► R 2 x R we define the maps ip x , p y , p z : SI x R — > R by 

!p(x,y,z) = (ip x (x,y, z),cp y (x,y, z), ip z (x,y, z)). 

We shall consider the collection T of all piecewise C° vector fields p : SI x R — » R 2 x R with the following 
property: there exists a finite number gi, . . . ,gk of functions in C 1 (Sl) such that the sets 



At 



{{x,y,z) : (x,y) € SI, gi(x,y) < z < g l+1 (x,y)} 



are nonempty and tp £ C 1 (yli,R 2 xR) for i = 0, . . . , k, where we put go — — oo and gk+i = +oo. 
Therefore, the discontinuity set of a vector field in T is contained in a finite number of regular surfaces. 

Let w G C 1 (rio) be a function such that J n \Vw\ 2 dxdy < +oo. The upper trace of w on S is 
denoted by w(x, 0+), and the lower trace by iy(a;, 0— ). Therefore, the approximate upper and lower 
limits w + (x,0) and u; - (ir,0) are given by 



w + (x, 0) = max{w(i, 0+), w(x, 0—)} 



and 



w (x,0) = min{w(x, 0+), w(x, 0— )}. 



A calibration for w is a bounded vector field tp E J 7 which is continuous on the graph of w and satisfies 
the following properties: 

(a) div(p — in the sense of distributions in SI x R ; 

(b) (ip x {x, y, z)) 2 + {tp y (x, y, z)) 2 < Aip z (x, y, z) at every continuity point (x,y,z) of tp; 

(c) (p x , p y )(x, y, w(x, y)) = 2Vjd(i, y) and p z (x, y, w(x, y)) = \Vw(x, y)\ 2 for every (x, y) € O ; 

(d) y^j ip x (x , y , z) dz^j + p y (x , y , z) dz^j < 1 for every (ar, y) G and for every ti,t% € R; 

(e) / ^(sC) 0, z)dz = and / P v (x, 0, z) rfz = 1 for every (x, 0) e S 1 . 

Jw-(x,0) Jw-( x ,0) 



The following theorem is proved in Q . 
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Theorem 2.1 If there exists a calibration ip for w , then w is a Dirichlet minimizer of the Mumford-Shah 
functional (1.1) in . 



If fl is a circle with centre on the x-axis, and w S C 1 (ilo) with J" n \S7w\ 2 dx dy < +00, then w 
satisfies the Euler conditions (a), (b), and (c) if and only if w has one of the following forms: 

J u(s,y) ify>0, 
l-u(x,y) + ci it y < 0, 



or 



/ \ \u(x,y) + c 2 ify>0, 
tw(ar,i/) = < (2.2) 
[it(aj,y) it y < 0, 

where u € C (Q) is harmonic with normal derivative vanishing on S and Ci , C2 are real constants. For 



our purposes, it is enough to consider the case ci = in (2.1) and C2 = 1 in (p.2j) 



3 A model case 



In this section we consider in (2.1) and in ( |2.2| ) the particular function u(x,y) = x and we deal with the 
minimality of the functions 

/ x / x if 3/ > 0, 
w(£>2/) : =< (3.1) 
[-a; if y < 0, 

and 

/ x jx + l ify>0, 
w(x,y):=<^ (3.2) 
I a; it y < 0. 

The aim of the study of these simpler cases (but we will see that they involve the main difficulties) is to 
clarify the ideas of the general construction. 



Theorem 3.1 Let w : R 2 — > R be the function defined by 

w(x,y) := 



x ify>0, 
x ify<0. 



Then every point (2:0,2/0) 7^ (0,0) has an open neighbourhood U such that w is a Dirichlet minimizer in 
U of the Mumford-Shah functional (1.1). 

Proof. The result follows by Theorem 4.1 of |jj if yo ^ 0. We consider now the case yo = , assuming 
for simplicity that xq > . We will construct a local calibration of w near (xq , 0) . Let us fix e > such 
that 

0< £ <^, 0< £ <i. (3.3) 
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y=-8 



y=5 



z=x 



Figure 1: Section of the sets A\, . . . ,A 5 at x — constant. 
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For < S < e we consider the open rectangle 

U := {(x,y) e R 2 : \x - x \ < e,\y\ < 6} 
and the following subsets of UxR (see Fig. 1) 



^3 

A 5 



= {(x,y, z) eU xR : x - a(y) < z < x + a(y)}, 

= {(x, y, z) e UxR : b + k(X) y < z <b + k(X) y + h} , 

= {(x, y,z) eU xR : —h < z < h}, 

= {(x, y, z) e UxR : -b + k(A) y-h<z<-b + k(A) y} , 

= {(x, y, z) e UxR : -x - a(-y) < z < -x + a(-y)}, 



where 



a 



(y) := v/4 £ 2 - ( £ - yf 



h := 



x — 3s 



K(A) := 4 - A' 



b :=2h + k(\)6, A:= 



l-4e 
2/i 



We will assume that 



5 < 



x — 3e 
8RA)|' 



(3.4) 



so that the sets A\, . . . , are pairwise disjoint. 

For every (x, y, z) £ UxR, let us define the vector ip(x, y, z) — (ip x , f v , f z )(x, y, z) € R 3 as follows: 



2(e - V) 



-2(z - x) 



^(e-yf + iz-x) 2 ' ^{e-yf + iz-xf 



0,A^ 



(/(»), 0,1) 
A 2 



0,A, 



2(e + y) 



2(z + z) 



V(e + y) 2 + (z + x) 2 ^/(e + y) 2 + (z + x) 



where 



(0,0,1) 



m := - j 



lj if (or, y,z) e A 1 , 

if {x,y,z) e A 2 , 
if {x,y,z) e A 3 , 
if (x,y,z) e At, 

r, 1 ^ if {x,y,z) e A 5 , 
otherwise, 



V t2 + (e - y) 2 



: (ft 



«(-y) 



v/i 2 + (e + yf 



:dt 
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Figure 2: Section of the set A\ at z = constant. 



Note that A\ U A 5 is an open neighbourhood of graph(w) n (f/xR). The purpose of the definition 
of ip in A\ and A 5 (see Fig. 2) is to provide a divergence free vector field satisfying condition (c) of 
Section 2 and such that 

<p v (x,0, z) > for \z\ < x, 
tp y (x,0, z) < for \z\ > x. 

These properties are crucial in order to obtain (d) and (e) simultaneously. 

The role of Ai and A4 is to give the main contribution to the integral in (e). To explain this fact, 
suppose, for a moment, that e — ; in this case we would have A\ = A 5 = and 

I <p y (x,0,z)dz = l, 

J — x 

so that the y-componcnt of equality (e) would be satisfied. 

The purpose of the definition of 93 in A3 is to correct the x -component of ip, in order to obtain (d). 

We shall prove that, for a suitable choice of 5, the vector field ip is a calibration for w in the rectangle 

U. 
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Note that for a given zgRwe have 

d x tp x {x,y, z) + d y tp y (x,y, z) = (3.5) 

for every (x,y) such that (x, y, z) G A± U A§ . This implies tp is divergence free in A\ U ^5 . Moreover 
div<^ = in the other sets Ai , and the normal component of tp is continuous across dAi : the choice of 
fc(A) ensures that this property holds for dAi and (see Fig. 3). Therefore tp is divergence free in 
the sense of distributions in UxH. 
On the graph of w we have 



tp(x,y,w(x,y)) 



( 2,0,1) if 2/ > 0, 
(-2,0,1) if 2/ < 0, 



so condition (c) is satisfied. 

Inequality (b) is clearly satisfied in all regions: the only non trivial case is A3 , where we have, using 

l/(v)l< 4( ° te)+ ° ( -» )) <^<2. 

xo — 3s xo — oe 

We now compute 

r x 

ip v (x,y, z) dz. (3-6) 



Let us fix y with \y\ < 5 . Since tp v {x, y, z) depends on z — x, we have 

/ <p»(x,y,z)dz = / ^(^y.x)^. (3.7) 

J x — a(y) J x 

Using ( |3.5[ ) and applying the divergence theorem to the curvilinear triangle 

T = {(£, V )eR 2 :H>x, r,<y, (e - t?) 2 + (x - £) 2 < 4e 2 } 

(see Fig. 4), we obtain 

+ " y, x) d£ = [ V tp x (x, n, x) d V = 2(y + s). (3.8) 



From (3.7) and (ET§), we get 



Similarly we can prove that 



f tpy(x,y,z)dz = 2(y + e). (3.9) 

J x-a(y) 



^ tpy(x,y,z)dz = 2(-y + e). (3.10) 
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Figure 4: The curvilinear triangle T . 



Using the definition of ip in A 2 , A 3 , A4 , we obtain 

/ <p v (x,y,z)dz = 1. 

J — X 

On the other hand, by the definition of / , we have immediately that 

I ip x (x,y lZ )dz = 0. 

J — X 



(3.11) 



(3.12) 



From these equalities it follows in particular that condition (e) is satisfied on the jump set S w fl U = 
{(x,y) eU :y = 0}. 

Let us begin now the proof of (d). Let us fix (x,y) G U . For every t\ < t 2 we set 

I(h,t 2 ) := / {<p x ,ip y ){x,y,z)dz. 
Jt 1 
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It is enough to consider the case —x — a(— y) < t\ < t 2 < x — a(y) . We can write 

I(h,t 2 ) = I(h,-x)+I(-x,x)+I(x,t 2 ), 
I(h,-x) = I(t 1 f\{-x + a{-y)),-x)+I{t 1 V {-x + a(-y)),-x + a(-y)), 
I(x,t 2 ) = I(x,h V (x - a(y))) + I(x - a(y),t 2 A (x - a(y))). 

Therefore 

I{hM) = J(-x, x) + 7(ti A {-x + a(-y)), -x) + I(x, t 2 V (x - a(y))) 

+ I(h V (-x + a{-y))M A (x - a{y))) - I(-x + a(-y),x - a(y)). (3.13) 

Let B be the ball of radius 4e centred at (0, — 4e) . We want to prove that 

I(x,t) e~B (3.14) 

for every t with x — a(y) < t < x + a(y). Let us denote the components of I(x,t) by a x and a v . 
Arguing as in the proof of (3.9), we get the identity 



a y = 2(e - y) - 2y/(t - x) 2 + (e - y) 2 < 0. (3.15) 

As | (/s 21 1 < 2 , we have also 

K) 2 < 4(t - x) 2 = (2(e - y) - a?) 2 - 4(e - y) 2 . 
From these estimates it follows that 

(a a: ) 2 + (a ! ' + 4£) 2 < 16e 2 , 



which proves ( 3.14 ). In the same way we can prove that 

I(t,-x)eB (3.16) 

for every t with —x — ot{—y) <t< —x + a(—y). 
If f(y) > 0, we define 

C := ([0, 2hf(y)] x[0, | - 2s}) U ({2hf(y)}x[Q, 1 - 4e]); 

if f(y) < 0, we simply replace [0,2hf(y)] by [2/i/(y),0]. From the definition of ip in A 2 , A$, A4, it 
follows that 

I(-x + a(-y),x - a(y)) = (2hf(y), 1 - 4e) (3.17) 

and 

%,s 2 )6C (3.18) 
for -a; + < si < s 2 < a; - a(y) . Let D := C - (2hf(y), 1 - 4e), i.e., 

£> = ([-2hf(y), 0] x [-1 + 4e, -\ + 2e]) U ({0}x [-1 + 4e, 0]), 
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for f(y) > 0; the interval \-2hf (y),0\ is replaced by [0, -2hf(y)] when f(y) < 0. From (|[I|), (|3.1l| ), 
( ^12|) , ( PI ), |U6| ), ( p^ ) and (|3.18|) we obtain 



J(ti, i 2 ) € (0,1) + 2B + D. 



As /(0) =0, we can choose 5 so that (^4) is satisfied and 

\2hf( y )\ = ^^\f(y)\<s 



(3.19) 



(3.20) 



for \y\ < S . It is then eas y to see that, by (3.2), the set (0, 1) + 2B + D is contained in the unit ball 
centred at (0,0). So that (|3.19|) implies (d). □ 



Remark. The assumption (xo,j/o) 7^ (0,0) in Theorem 3.1 cannot be dropped. Indeed, there is no 
neighbourhood U of (0,0) such that w is a Dirichlet minimizer of the Mumford-Shah functional in U . 

To see this fact, let tp be a function defined on the square Q = (— 1, 1) x (—1, 1) satisfying the boundary 
condition ip — w on dQ and such that = ((—1, —1/2) U (1/2, l))x{0}. For every e, let ip e be the 
function defined on Q £ = eQ by ip e (x,y) := sip(x/e, y/e) . Note that ip s satisfies the boundary condition 
tj} e = w on dQ e . Let us compute the Mumford-Shah functional for if> e on Q e : 

[ \Vip e \ 2 dxdy + H 1 (S i , e ) = e 2 [ \ \7^\ 2 dx dy + e. 

I iVwfdxdy + H 1 ^) = 4e 2 + 2e, 
JQe 

[ iVipefdxdy + HXS^jK [ \Vw\ 2 dxdy + H 1 {S W ) 



Since 



we have 



for e sufficiently small. 



The construction shown in the proof of Theorem 3.1 can be easily adapted to define a calibration for 
the function w in (3.2). 



Theorem 3.2 Let 



R 2 — ► R be the function defined by 
w(x,y) :-- 



1 ify>0, 
ify < 0. 



Then every point (xo,yo) € R- 2 has an open neighbourhood U such that w is a Dirichlet minimizer in 
U of the Mumford-Shah functional (1.1). 

Proof. The result follows by Theorem 4.1 of Q if yg ^ . We consider now the case yo — 0; we will 
construct a local calibration of w near (xq, 0) , using the same technique as in Theorem 3.1 We give only 
the new definitions of the sets A\ , . . . , A5 and of the function ip , and leave to the reader the verification 
of the fact that this function is a calibration for suitable values of the involved parameters. 
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Let us fix e > such that 



° <£< 24' 



0<£< h- 



For < 5 < e we consider the open rectangle 

U := {(x,y) g R 2 : \x - x \ < e, \y\ < 5} 
and the following subsets of (JxR 



A 2 
As 
A 4 
A 5 



{(x, y,z) eU xR : x + 1 - a(y) < z < x + 1 + a(y)}, 
{(x, y, z) g C/xR : & + k(A) y + 3ft < z < b + k(X) y + 4ft}, 
{(a;, j/, z) g C/xR : .t + 3e + 2ft < z < x + 3e + 3ft}, 
{(x, y, z) g UxH : b + k(\) y < z <b + k(X) y + ft}, 
{(x, y,z) g £7 xR : x - a(-y) < z < x + a(-y)}, 



where 



a 



(y) := ^Ae 2 ~ (e-y) 2 , 



ft := 



1 - 6e 



/x\ Al , 1 — 4e 

k(A) :=---, 6 := x + 3e + k(A) d, A := — — . 
4 A 2/z 



We will assume that 



6 < 



l-6s 

ToRa)|' 



so that the sets A\, . . . ,A$ are pairwise disjoint. 

For every (x, y, z) G [7xR, let us define the vector ip(x, y, z) g R 3 as follows: 



2(e-y) 



-2(z-z- 1) 



V( £ -2/) 2 + (z-x-l)2' ^( £ -yr + (z-x-l)2 



0,A,£ 



(/(y),o,i) 

0,A,^ 



2(e + y) 



2(z-x) 



^TWTi^xY' v / (£ + 2/) 2 + (^-x) 2 ' 



(0,0,1) 



,lj if (jc, y, z) g A x 

if (x,y,z) g A 2 
if (x,y,z) g A 3 
if (x,y, z) g A 4 

if (x,y, z) g A 5 
otherwise, 
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where 




dt+ / = =dt 



Vt 2 + (e ~ V? Jo Vt 2 +¥+V) 



for every \y\ < 5. □ 



4 The general case 



In this section we denote by O a ball in R 2 centred at (0,0) and we consider as u in (2.1) and in (2.2) 
a generic harmonic function with normal derivative vanishing on S . We add the technical assumption 
that the first and second order tangential derivatives of u are not zero on S . 

Theorem 4.1 Let u : il — > R be a harmonic function such that d y u(x,0) — for (x,0) £ SI, and let 
w : SI — > R be the function defined by 



w(x,y) 



u(x, y) for y > 0, 
-u(x, y) for y < 0. 



Assume that uq := u(0,0) 7^ 0, d x u(0,0) 7^ 0, and <9 2 u(0,0) 7^ 0. Then there exists an open neighbour- 
hood U of (0, 0) such that w is a Dirichlet minimizer in U of the Mumford-Shah functional fil.^ ). 

Proof. We may assume u(0, 0) > and d x u(0, 0) > 0. We shall give the proof only for <9 2 u(0, 0) > 0, 
and we shall explain at the end the modification needed for <9 2 u(0, 0) < 0. Let v : SI — > R be the 
harmonic conjugate of u that vanishes on y = 0, i.e., the function satisfying d x v(x,y) — — d y u(x, y), 
d y v{x, y) = d x u{x, y) , and v{x, 0) = 0. 

Consider a small neighbourhood U of (0,0) such that the map $(x, y) := (u(x,y),v(x,y)) is invertible 
on U and d x u > on U. We call ^ the inverse function (u,v) 1— > (£(u, v), 77(11, v)) , which is defined in 
the neighbourhood V := $(£/) of (uq,0)- Note that, if U is small enough, then ri(u,v) = if and only 
if v = 0. Moreover, 

D* = ( ) = * ( f U f V ) , (4.1) 

V 9 u ri d v rj J |Vu| 2 \ a y u d y v ) 

where, in the last formula, all functions are computed at (x, y) = ^(u, v) , and so <9 U £ = d v r], d v £ = —dui] 
and d u r](u, 0) = , d v rj(u, 0) > . In particular, £ and r\ are harmonic, and 

^T7(tt,0) = 0, d^(u,0) = 0. (4.2) 



On U we will use the coordinate system (u,v) given by $. By (4.1) the canonical basis of the tangent 
space to U at a point (x, y) is given by 

VU VV (4.3) 



For every (u,v) G V , let G(u,v) be the matrix associated with the first fundamental form of U in the 
coordinate system (u,v), and let g(u,v) be its determinant. By (4.1) and (4.3), 

a=((9uvf + (d v v) 2 ) 2 = W ^ W - (4.4) 
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We set j(u, v) := f/g(u,v). 

The calibration (p(x,y, z) on (7xR will be written as 



l 



(u(x,y),v(x,y),z). 



l 2 (u(x,y),v(x,y)) 
We will adopt the following representation for <j> : VxR — ► R 3 : 

w, 2) = w, z)t„ + 4> v {u, v, z)t v + 4> z (u, v, z)e z , 



(4.5) 



(4.6) 



where e z is the third vector of the canonical basis of R 3 , and r„ , t v are computed at the point ^(u, v) . 
We now reformulate the conditions of Section 2 in this new coordinate system. It is known from Differ- 
ential Geometry (see, e.g., |l], Proposition 3.5]) that, if X = X w t u + X v t v is a vector field on U , then 
the divergence of X is given by 



divX = -1(9 U ( 7 2 X«) + d v ( 7 2 X v )). 
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(4.7) 



Using (4.3), (4.4), (fa ,5|), (l.Q, and (1.7) it turns out that 9? is a calibration if the following conditions 
are satisfied: 

(a) d u (j) u + d v <j? + d z (j) z = for every (u, v, z) G VxR; 

(b) (c/) u (u,v,z)) 2 + (c/) v (u,v,z)) 2 < 44> z (u,v,z) for every (u, v, z) G VxR; 

(c) (j) u (u,v,±u) = ±2, <t>°(u, v, ±u) = 0, and ^> z (u,w,±u) = 1 for every (u,v) £ V; 



(d) 



b u (u,v, z) dz + 



'1 



(u,v, z) dz \ < j 2 (u, v) for every (u,v) G V , ii , £2 € R; 



(e) / 4> u (u,ti,z)dz = and / v (u,O, z) = 7(w,0) for every (u, 0) e V. 

J —u J —u 

Given suitable parameters £ > 0, h > 0, A>0, that will be chosen later, and assuming 

V = {{u,v) : \u-u \ < S,\v\ < 5}, 
with 5 < e, we consider the following subsets of FxR 



(4.8) 



Ax 
A 2 

A 4 
A 5 



= {{u, v, z) £ VxR : u — a(v) < z < u + a(v)}, 

= {(u,v,z) G VxR : 3h + 0(u,v) < z < 3h + (3(u,v) + 1/A} , 

= {(u, v, z) G VxR : — h< z < h}, 

= {(u,v, z) G VxR : -3h + 0(u,v) -1/X<z< -3h + 0(u,v)} , 

= {(u, v , z) £ V xR : —u — a(— v) < z < — u + a(—v)}, 



where 



y(v) := \/4e 2 - (e - v) 2 



and P is a suitable smooth function satisfying j3(u, 0) = 0, which will be defined later. It is easy to see 
that, if e and h are sufficiently small, while A is sufficiently large, then the sets A%, . . . , A5 are pairwise 



14 



disjoint, provided 5 is small enough. Moreover, since 7(1*, 0) = d v r)(u, 0) > 0, by continuity we may 
assume that 



•y(u, u) > 128e and d v ri(u,v) > 8e 

for every (u, v) € V ■ 

For (u, w) e V and z e R the vector (f)(u, v, z) introduced in (|4.5| ) is defined as follows: 



2(e - w) 



2(z - it) 



\J{e — v) 2 + (z — u) 2 tJJe — v) 2 + (z — u)' 



— Acr(u, v) 

f(v)r u + e 

— \<j(u, V 



\/{u- a) 2 + v 



-.t u + Aer(w, v) 



^J(u — a) 2 + v 2 



y/ (u — a) 2 + v 2 
2(e + w) 



t u + Xa(u, v) 



\J{u — a) 2 + v 2 
2{z + u) 



y/(e + v) 2 + (z + u) 2 ^/(e + v) 2 + (z + u) 2 



in A x , 

r v +iie z in A 2 , 

in A 3 , 
t v + fie z in At, 

e z in A 5) 

otherwise, 



where 



a < uq — 11(5, > 



/(«) 



1 

h 



a(v) 



(e-v) 



^t 2 + (e- vf 



:dt 



a(—v) 



(e + v) 



yjt 2 + (e + vf 



: (It 



a(u, v) := ±7(0 + ^/(w - a) 2 + v 2 , 0) - 2e. 
We choose /3 as the solution of the Cauchy problem 

\a(u, v)(—v d u j3 + (u — a)d v (3) = (/x — l)y/(u — a) 2 + v 2 , 

< 

J(u,0)=0. 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



Since the line v = is not characteristic for the equation near (uo,0), there exists a unique solution 
/3 G C°°{V) , provided V is small enough. 

In the coordinate system (u,v) the defi nition of the field <p in Ax, A3, and A5 is the same as the 
definition of ip in the proof of Theorem 3.1. The crucial difference is in the definition on the sets A 2 and 
A4 , where now we are forced to introduce two new parameters a and /1 . Note that the definition given 
in Theorem 3.1 can be regarded as the limiting case as a tends to +00. 
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By direct computations it is easy to see that <fi satisfies condition (a) on A\ and A$ . Similarly, the 
vector field 

/ v u — a \ 



\ yj(u — a) 2 + v 2 yj {u — a) 2 + v 2 I 

is divergence free; since (u — a) 2 + v 2 is constant along the integral curves of this field, by construction 
the same property holds for a , so that </> satisfies condition (a) in A2 and A4 . 
In A3 , condition (a) is trivially satisfied. 

Note that the normal component of <f> is continuous across each dAi : for the region A3 this continuity 
is guaranteed by our choice of /3. This implies that (a) is satisfied in the sense of distributions on V^xR. 
In order to satisfy condition (b), it is enough to take the parameter \i such that 



A 5 



C7 2 (u, v) < [I 



for every (u, v) € V , and require that 



\f(v)\ < 2- 



(4.13) 



Since 



*(v) + a(-v) As 
h h 



(4.14) 



inequality ( 4.13| ) is true if we impose 



2e< h. 



Looking at the definition of <f> on A\ and A^ , one can check that condition (c) is satisfied. 
Arguing as m the proof of ( |3~7o| ), fl3~12| ) in Theorem pl|, we find that for every (u, v) G V 



— u+a{ — v) r-h 

(j) w (u,v, z) dz + J (j) w (u,v, z) dz + j <j) u (u,v, z) dz = 0, 



-u+a(— v) ph c-u 

<fi v (u,v, z) dz + / <fi"(u, v, z) dz + / cf> v (u,v, z) dz = 4e. 

■j, J —h Ju — a(v) 



Now, it is easy to see that 



0"(w, v, z) dz = — 2a(u, v) 



^{u-a) 2 + 



(4.15) 



since for v — we have 



condition (e) is satisfied. 



4> v (u, v, z) dz = 4e + 2a(u, v) 



y/(u — a) 2 + V 2 



a(u,0) = - 7 («,0)-2e, 



(4.16) 
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By continuity, if 5 is small enough, we have 

f u 7 
J v {u,v,z)dz> -7(«,») (4.17) 

for every (u, v) G V . 

From now on, we regard the pair (0", (ff) as a vector in R 2 . To prove condition (d) we set 

I eja (u,v,s,t) := J ((/) u ,(j) v )(u,v,z)dz 

for every (u,v) G V, and for every s,t G R. We want to compare the behaviour of the functions |/ e , | 
and 7 2 ; to this aim, we define the function 

d eta (u, v, s, t) := \I s , a (u, v, s, i)| 2 - -y 2 (u, v). 

We have already shown (condition (e)) that 

d e , a (u,0,-u,u) = 0. (4.18) 

We start by proving that, if V is sufficiently small, condition (d) holds for every (u,v) € V , for t\ close 
to — u and t2 close to u. Using the definition of cj)(u, v, z) on A\ and A5, one can compute explicitly 
d £ ^ a (u, v, s, t) for \s + u\ < a(— u) and for |t — m| < a(v) . By direct computations one obtains 



V Ol ,,td Bl a(«,0,-«,«) = (4.19) 



for (tt,0) G 



We now want to compute the hessian matrix V 2 s t d £ ^ a at the point (mq, 0, — ito, wq) ■ By (4.11) and 



(4.4), after some easy computations, we get 

d%<r(u,0) = — - — -d u j{u,0) = — - — -d v d u r){u,0). 
2{u — a) 2(u — a) 

Using this equality and the explicit expression of d e ^ a near {uq, 0, — uq, uq) , we obtain 

5 2 d Eia (uo,0,-u ,uo) = -7 rz(d v r)(u ,0) - 4e) + 

(u - ay 

2 

-d v T](u o ,0)d v d u ri(u o ,0) - d 2 {j 2 )(u ,0). 



Uq — a 



Since rj and 7 do not depend on a and e, for every e satisfying (4.9) we can find a so close to uq that 

a 2 d £ , a (w o ,0,-M O ,uo) < 0. (4.20) 

Moreover, we easily obtain that 

d 2 d E:a (uo,0,~u O7 uo) = d 2 d £ , a (u o ,0,-u o ,uo) = 8- -d v r](u , 0), 

e 

4 

d v d t d e , a (u o ,0, -u ,u ) = d v d s d e} a(u o ,0, -u ,u ) = (d v r](u o ,0) - 4e), 

uo — a 

dtd s d E . a {u o ,0, -u 0) u ) =8. 
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By the above expressions, it follows that 

d%d eta d v d t d E! a 



det 
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(u ,0, -Uo,Uq) 



(u - a) : 



d v rj(u o ,0)(d v r](uo,0) - 4e) + 



ci(g) 
ito — a 



c 2 (e), 



where c\(e) , 02(e) are two constants depending only on e. Then, if e satisfies ( |4.9| ), a can be chosen so 
close to uq that 



det 



(u ,0, -u ,u ) > 0. 



dyd Eta d v d t d Sia 
d v d t d e , a d?d e , a 

At last, the determinant of the hessian matrix of d £<a at (uq, 0, — Uq, uq) is given by 



(4.21) 



detV L,t4,a("0,0,-U ,U ) = 



1 



32 



Uq — a 



-(d v r)(u , 0)yd v d u -q{u ,0)(d v T](u , 0) - 4e)— + c 3 (e) 



where 03(e) is a constant depending only on e. Since, by (4.1), 



d v d u rj(u o ,0) 



(d x u(o,o)) 3 ' 

given e satisfying ( |4.9| ), we can choose a so close to uq that 

det V^ Sjt d St a(u , 0, -u , u ) < 0. 



(4.22) 



By ( 4.20|) , ( 4.21 ), and ( 4.22j ), we can conclude that, by a suitable choice of the parameters, t he he ssian 
matrix of d £J1 (with respect to v,s,t) at (uq, 0, — uq, uo) is negative definite. This fact, with (4.18) and 
(4.19), allows us to state the existence of a constant r > such that 



d e ,a(u, V, S,t) < 



(4.23) 



for \s + u \ < t, \t— u \ < t, (u,v) £ V , v 7^ 0, provided V is sufficiently small. So, condition (d) is 
satisfied for \ti + u \ < r and |*2 — «o| < T - We can assume 6 < t < a(v) for every (u,v) S V. 

From now on, since at this point the parameters e , a have been fixed, we simply write / instead of 
Ie,a ■ We now study the more general case \t\ + u\ < a(—v) and \t<2 — u\ < a(y) . 

Let us set 

mi(u, v) := max{|/(it, v, s, t)\ : \s + u\ < a(—v), \t — u\ < a(v), \t — uq\ > r} . 

By the definition of Ax,... , for p = a (5) + S we have (0", (jf) = on (Vx [uq — p,u + p\) \ A\ and 
(Vx [— u — p, —u + p]) \ A 5 . This implies that 

mx(u, v) :— max{|J(u, v, s,t)\ : \s + u \ < p, r < \t — u \ < p} 

for (u,v) G V . The function mi, as supremum of a family of continuous functions, is lower semicon- 
tinuous. Moreover, mi is also upper semicontinuous; indeed, suppose, by contradiction, that there exist 
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two sequences (u n ), (v n ) converging respectively to u, v, such that (mi(%, v n )) converges to a limit 
I > mi(u,v); then, there exist (s n ) , (t n ) such that 



|s„ + u n | < a(-v n ), \t„ - u n \ < a{v n ), \t n - u \ > r, 



(4.24) 



and mi(it n ,t> n ) = tm, s n , t n )\ . Up to subsequences, we can assume that (s„), (t n ) converge 

respectively to s, t such that, by ( 4.24 ), 



|s + u\ < a(— v), \t — u\ < a(v), \t — uq\ > r; 



hence, we have that 



mi(u,v) > \I(u,v, s,t)\ = lim \I(u n ,v n ,s n ,t n )\ — I > mi(u,v), 

n — >oc 

which is impossible. Therefore, mi is continuous. 

Let B be the open ball of radius 4e centred at (0, — 4e) . Arguing as in (3.14), we can prove that 

I{u,v,u,t)eB (4.25) 

whenever < \t — u\ < a(v) . In the same way we can prove that 

I[u, v, s, —it) £ B (4.26) 

for < \s + u\ < a(—v) . We can write 

I(u, v, s, t) = I(u, v, s, —u) + I(u, v, —u, u) + I(u, v, u, t). (4.27) 



So, for \s + u\ < a(-v), \t-u\ < a(v), and \t - u \ > r, by ( ^.26| ), ( |4.15| ), ( [4.16[ ), and ( ^.25| ), we obtain 
that 

I(u, 0, s, t) € (0, 7(u, 0)) + B + B, 



hence, by (4.9), I(u, 0, s, t) belongs to the open ball of radius j(u, 0) centred at (0, 0) , and so, mi (it, 0) < 
7(m, 0) . By continuity, if V is small enough, 



mi(u, v) < j(u, v) 



(4.28) 



for every (u, v) € V . 

Analogously, we define 



7712(11, v) :— max{|/(it, v, s, t)\ : \s + u\ < a{~v), \s + uq\ > r, \t — u\ < a{v), } . 
Arguing as in the case of mi , we can prove that, if V is small enough, 

m 2 (u,v) < j(u,v) (4.29) 

for every (w, v) € V . 

By ( 1.28| ), (4.29), and (4.23), we can conclude that I{u,v,t\,t2) belongs to the ball centred at (0,0) 
with radius 7(77,7;), for \t\ + u\ < a(—v) and \t 2 — u\ < a(v) . More precisely, let E(u,v) be the 
intersection of this ball with the up per h al f plan e boun ded by the horizontal straight line passing through 
the point (0, §7(77,7;)): by J4.27p , ( [4.17D , ( |4.25| ), (f4.26| ), and (Q, we deduce that 

I(u,v,t 1 ,t 2 ) S E(u,v) (4.30) 
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for \ti + u\ < a(—v) and |ta — u\ < a(v) . 

We can now conclude the proof of (d). It is enough to consider the case — u — a(—v) < t\ < t 2 < 
u + a(v) . We can write 



I(u,v,t 1 ,t 2 ) = I(u, v,ti A (— u + a(— v)),t% V (u — a(v))) + 
+I(u, v, ti V (-u + a(-v)), t 2 A (u - a(v))) - 
—I(u, v, —u + a{— v), u — oi(v)). 



By ( |4.30| ), it follows that 

I(u, v, ti A (—u + a(~v)),t 2 V (it - G S(tt, u). 

Let Ci(u,f) be the parallelogram having three consecutive vertices at the points 

(— v, u — a) 



(2hf(v),0), (0,0), a(u,v) 

let C2 (u, v) be the segment with endpoints 

(2ft/(«),0), (2ft/(«),0) + 2CT(«,«) 

and let C(w, u) := Ci(u, w) U C2(m, w) . 

From the definition of tp in A 2 , A3, A4, it follows that 



y/(u — a) 2 + v 2 

(—v, u — a) 
yj(u — a) 2 + v 2 



I(u, v, —u + a(—v),u — a>(v)) — (2hf(v), 0) + 2a(u, v) 



(—v, u — a) 
y/(u — a) 2 + V 2 



and 



I(u, v, si, s 2 ) € C(u, v) 
for —u + a{— v) < Si < S2 < u — . Let 

L>(u,w) := C(«,w) - (2/i/(w),0) - 2cr(u,w)- 



(— v , u — a) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



y/(u — a) 2 + V 2 

From ( |4.31[ ), ( |4.32[ ), ( |4.33| ), and ( |4.34D we obtain 

7(u, w, £1, £ 2 ) e u) + D(u, u). (4.35) 

As |«| < S < 105 < u — a by ( |4. 1C| ) , the angle that the segment C 2 (u,v) forms with the vertical is less 
than arctan(l/10) . Moreover, we may assume that the lenght 2a(u,v) of the segment C 2 (u,v) is less 
t han y(u, v) ; in deed, this is true for v = and, by continuity, it remains true if 6 is small enough. By 
(4.9) and ( 4.14 ), we have also that \2hf(v)\ < j(u,v)/16. Using these properties and simple geometric 
considerations, it is possible to prove that E(u,v) + D(u,v) is contained in the ball with centre (0,0) 
and radius j(u,v). This concludes the proof of (d). 
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If d 2 u(0, 0) < 0, it is enough to change the definition of cj) in the sets A2 and A4 , as follows: 

v . , a — u 



\a(u, v) 



where a > uq + 118 and 



□ 



i/(a — u) 2 + v' 



-.t u + Xa(u,v) 



y/(a — u) 2 + v 2 



Theorem 4.2 Let w : — > R 6e a harmonic function such that d y u(x,0) — for (x,0) G f2, and iei 
u; : 51 — > R 6e i/ie function defined by 



)(x,y) := 



u(x,y) + l fory>0, 
u(x,y) fory<0. 



Assume that d x u(0,0) 7^ and d 2 u(0,0) 7^ 0. Then there exists an open neighbourhood U of (0,0) such 



that w is a Dirichlet minimizer in U of the Mumford- Shah functional (1.1) 



Proof. We will write the calibration (p as in (4.5) and we will adopt the representation (4.6) for 
4>. We will use the same technique as in Theorem 4.1. We give only the new definitions of the sets 
Ax,... ,A$ and of the function (j> when d x u(0, 0) > and d 2 u(Q, 0) > 0, and leave to the reader the 
verification of the fact that this function is a calibration for suitable values of the involved parameters. 



The case d 2 u(0, 0) < can be treated by the changes introduced at the end of Theorem 4.1 
Let uq := u(Q, 0) . Given e > 0, h > 0, A>0, and assuming 

V := {(u,v) : \u-uq\ < S,\v\ < S}, 

we consider the following subsets of VxR 



Ax 


= {(«, 


v, 


z) e VxR 


U + 1 — a(v) < z < u + 1 + ce(v)}, 


A 2 


= {(«, 


v, 


z) e VxR 


5h + (3(u, v) < z <5h + 0(u, v) + 1/A}, 


A 3 


= {(«, 


v, 


z ) e VxR 


2h < z < Ah}, 


A A 


= {(«, 


v, 


z) e VxR 


h + f3(u,v) < z < h + j3(u,v) + 1/A}, 


A 5 


= {(«, 


v, 


z) G VxR 


u — a(— v) < z < u + a(— u)}, 



where 



a(i>) := \J Ae 2 - (e - w) 2 
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and (3 is a suitable smooth function satisfying /3(u, 0) = , which will be defined later. For (u, v) S V 
and zeR the vector <fi(u, v, z) is defined as follows: 



2(e-v) 



2(z-u- 1) 



-T v + e z in A u 



y/(e - v) 2 + (z - u - l) 2 y/(e - v) 2 + (z-u - 1) : 

V . , . u — a . , 

— Act(m, i>) ;/ Nn — - t u + Act(m, u)— — xn o r t ,+/xe z m A 2 , 



\JJu — a) 2 + v 2 



y/(u-a) 2 +v 2 



f(v)r u + e 
— Acr(u, v) 



in A 3 , 



y/ {u — a) 2 + V 2 
2{e + v) 



t u + Xa(u, v) 



yf{u-a) 2 +v 2 
2(z - u) 



t v + [ie z in A4, 



^(e + v) 2 + (z-u) 2 y/(e + v) 2 + (z - uj- 



--r v + e z 



in As, 



otherwise, 



where a < uq — 11(5, [l > 0, 



m ■= t 



Vt a + (e-«) : 



V* 2 + (e + vf 



a(u, v) := \-y{a + ^/(w - a) 2 + v 2 , 0) - 2e, 



and /3 is the solution of the Cauchy problem (4.12). □ 
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